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Abstract. This contribution discusses isocurvature modes, in particular the non-Gaussianities 
of local type generated by these modes. Since the isocurvature transfer functions differ from the 
adiabatic one, the coexistence of a primordial isocurvature mode with the usual adiabatic mode 
leads to a rich structure of the angular bispectrum, which can be decomposed into six elementary 
bispectra. Future analysis of the CMB data will enable to measure their relative weights, 
or at least constrain them. Non-Gaussianity thus provides a new window on isocurvature 
modes. This is particularly relevant for some scenarios, such as those presented here, which 
generate isocurvature modes whose contribution in the power spectrum is suppressed, as required 
by present data, but whose contribution in the non-Gaussianities could be dominant and 
measurable. 



1. Introduction 

Inflation is currently the best candidate to explain the generation of primordial perturbations 
(see e.g. [1] for a recent pedagogical introduction), but current observations cannot point to 
a specific scenario. The hope is thus that future data will enable us to discriminate between 
various models. In this respect, an important distinction is between single-field and multiple- 
field models. The detection of even a tiny fraction of isocurvature mode in the cosmological 
data would rule out single- field inflation, which predicts only adiabatic perturbations. Another 
signature that could distinguish multiple-field models from single-field models is a detectable 
primordial non-Gaussianity of the local type. Similarly to isocurvature modes, a detection of 
local primordial non-Gaussianity would rule out all inflation models based on a single scalar 
field, since they generate only unobservably small local non-Gaussianities. 

This contribution, based on the recent works [21 El SI E] , discusses how isocurvature modes 
could affect non-Gaussianities and their specific signature in the Cosmic Microwave Background 
(CMB) fluctuations. In particular, it is shown how the bispectrum generated by the adiabatic 
mode together with one isocurvature mode leads to a total angular bispectrum which can be 
decomposed into six distinct components: the usual purely adiabatic bispectrum, a purely 
isocurvature bispectrum, and four other bispectra that arise from the possible correlations 
between the adiabatic and isocurvature mode. Because these six bispectra have different shapes 
in angular space, their amplitude can in principle be measured in the CMB data. 

This analysis opens a new window on isocurvature modes, especially important to test models 
where the isocurvature contribution is suppressed in the power spectrum but not in the non- 
Gaussianities. Examples of such models are presented at the end of this contribution. 



2. Iso curvature modes 

At the time of last scattering, the main components in the Universe are the CDM (c), the 
baryons (b), the photons (7) and the neutrinos (y). All these components are characterized 
by their individual energy density contrasts 5{. The most common type of perturbation is the 
adiabatic mode, characterized by the condition 
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$c = $b = jd v = -5 y , (1) 

which means that the number of photons (or neutrinos, or CDM particles) per baryon does not 
fluctuate. Assuming adiabatic initial conditions is natural if all particles have been created by 
the decay of a single degree of freedom, such as a single inflaton, and, so far, the CMB data are 
fully compatible with purely adiabatic perturbations. 

However, other types of perturbations can be included in a more general framework. In 
addition to the adiabatic mode, one can consider four distinct isocurvature modes [6]: the CDM 
isocurvature mode, the baryon isocurvature mode, the neutrino density isocurvature mode and 
the neutrino velocity isocurvature mode. The first three isocurvature modes are characterized 
by 

s x = r^—tx X = {c, b, vd] (2) 
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with wx = Px/px- As for the neutrino velocity isocurvature mode, it is characterized by a non 
vanishing "initial velocity", compensated by the velocity of the photon-baryon plasma so that 
the total momentum density is cancelled, while the energy densities satisfy ([I]). 

In the following, these five modes will be denoted collectively as X . In the context of 
inflation, a necessary, although not sufficient, condition for at least one of these isocurvature 
modes to be produced is that several light degrees of freedom exist during inflation. Moreover, 
since the adiabatic and isocurvature modes can be related in various ways to these degrees of 
freedom during inflation, one can envisage the existence of correlations between these modes |7J. 

The various modes lead to different predictions for the CMB temperature and polarization. 
Let us consider for instance the temperature anisotropics, which can be decomposed into 
spherical harmonics: 
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At linear order, the multipole coefficients ai m are related to a general primordial perturbation 
consisting of the superposition of several modes, via the expression 
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where g((k) denotes the transfer function associated with the mode X 1 . As a result, the total 
angular power spectrum is given by 
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where the primordial power spectra Pij{k) are defined by 

(X I (k 1 )X J (k 2 )} = (2vr) 3 5(ki + k 2 )P /J (fci) . (6) 



Since the various transfer functions are different, this leads to different predictions for the CMB 
angular power spectrum. This is illustrated in Fig. [T] (left panel), where the angular power 
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Figure 1. Angular power spectra (multiplied by 1(1 + 1)) for the temperature (left) and 
polarization (right) obtained from purely adiabatic or purely isocurvature initial conditions. 
The amplitude and spectral index of the primordial power spectrum, as well as the cosmological 
parameters, on which the transfer functions depend, correspond to the WMAP7-only best-fit 
parameters. 



spectra produced separately by the various modes are plotted, assuming the same primordial 
power spectrum. The only exception are the CDM and baryon isocurvature modes which give 
exactly the same pattern, up to the rescaling Sb = (flb/Qc) S c where £lb and £l c denote, as usual, 
the present energy density fractions, respectively for baryons and CDM. 

We infer from CMB observations that the "primordial" perturbation is mainly of the adiabatic 
type. However, this does not preclude the presence, in addition to the adiabatic mode, of 
an isocurvature component, with a smaller amplitude. Precise measurement of the CMB 
fluctuations could lead to a detection of such an extra component, or at least put constraints on 
its amplitude. For example, constraints on the CDM isocurvature to adiabatic ratio, 

« = (7) 

based on the WMAP7+BAO+SN data, have been published for the uncorrelated and fully 
correlated cases (the impact of isocurvature perturbations on the observable power spectrum 
indeed depends on the correlation between adiabatic and isocurvature perturbations, as 
illustrated in [8]). In terms of the parameter a = a/(l + a), the limits given in [9] are 

a < 0.064 (95%CL), ax < 0.0037 (95%CL) , (8) 

respectively for the uncorrelated case and for the fully correlated case. 

3. Generalized angular bispectra 

Let us now turn to the CMB non-Gaussianities. The angular bispectrum corresponds to the 
three-point function of the multipole coefficients: 

hhh = \ a hm 1 ai2m2 a hm 3 ) ■ {&) 

Substituting the expression Q into the angular bispectrum, one finds 

K?z 2m3 = Qwi?r A K hh , (io) 

where the first, purely geometrical, factor is the Gaunt integral 

^™ m3 = / ^ W*) W*) W*0. (H) 



while the second factor, usually called the reduced bispectrum, 
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depends on the bispectra of the primordial X 1 : 

(X I (k 1 )X J (k 2 )X K (k 3 )) = (27r) 3 5(^ l k i )B IJK (k 1 , k 2 , k 3 ) . (13) 

We now need to specify the primordial bispectra B IJK . For purely adiabatic perturbations, 
the local bispectrum is expressed as the square of the power spectrum (symmetrized over k\ , k 2 
and k 3 ). Here, we consider the generalization 

B IJK (h, k 2 , ks) = f^ K Pdk 2 )Pdh) + fif'PdhWks) + f£i IJ P c (h)Pdh) , (14) 

where the coefficients f^ K satisfy the condition 

f,JK _ f,KJ n ,x 

The above expression is the natural outcome of a generic model of multiple-field inflation. Indeed, 
allowing for several light degrees of freedom during inflation, one can relate, in a very generic 
way, the "primordial" perturbations X 1 (defined during the standard radiation era) to the 
fluctuations of light primordial fields (fi a , generated at Hubble crossing during inflation, so that 
one can write, up to second order, 

X 1 = Ni 5^ + l -Ni b 54> a 5ct> h + ... (16) 



where the 5(p a can usually be treated as independent quasi-Gaussian fluctuations, i.e. 

<^(k)5/(k')) = (27rf5 ab P 5<l> (k)5(k + -k'), P S4> (k) = 27r 2 £T 3 (^] , (17) 



where a star denotes Hubble crossing time. The relation (16) is very general, and all the details 
of the inflationary model are embodied by the coefficients and NL. Substituting (16) into 
( 13 ) and using Wick's theorem, one finds that the bispectra Bjjk can be expressed in the form 

B IJK (k l7 k 2 , k 3 ) = X^Ps^Ps^h) + \ J ' KI P S< p(ki)P d<i> (k 3 ) + \ K > IJ ' P S4> {k x )P s<i) {k 2 ) , (18) 

with the coefficients 

X I,JK = § ac S bd N I bN J N K (lg) 

(the summation over scalar field indices a, b, c and d is implicit), which are symmetric under 
the interchange of the last two indices, by construction. Since the adiabatic power spectrum is 
given by 

P C = (5 ab N< N ( b )P s ^ = AP S(f> , (20) 



one obtains finally (14) with 

/n'l* = A^/M 2 , (21) 

where it is implicitly assumed that the coefficients are weakly time dependent so that the 
scale dependence of A 2 can be neglected. 



After substitution of (14) into (12), the reduced bispectrum can finally be written as 
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where each contribution is of the forrrQ 

b hH = 3 /„ r*dral h (rW 2 (r)f3«(r), (23) 

with 
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k 2 dkjt(kr) gf(k), P{(r) = - / k 2 dkjt(kr) gj(k) P c (k) . (24) 
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4. Observational prospects 

For simplicity, we assume that the primordial perturbation is the combination of the dominant 
adiabatic mode with a single isocurvature mode. In this case, the total bispectrum is 
characterized by six parameters, which we now denote /W, 

K hl3 = F« + 2fM tffi + e + f*« bf;fi + b fff is + f s,ss 6 sg 

= E/ ( HV ( 25 ) 
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where the index i varies between 1 to 6, following the order indicated in the upper line. Note that, 

because of the factor 2 in front of f^^ s and f s, ^ s , we define b^ 2 j = ^ili^h and ; j = ^b^fofs 
whereas there is no such factor 2 for the other terms. 

To estimate these six parameters, given some data set, the usual procedure is to minimize 

(B obs - f {l) B^), (B obs - f {i) B^)\ , (26) 
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where B is a short notation for the angle- averaged bispectrum 
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For an ideal experiment (no noise and no effects due to the beam size) without polarization, the 



scalar product in ( 26 ) is defined by 



(B,B')= ]T Bhhl f' 1 ^ . (28) 
h<h<h a hhh 

with the variance 

<hh = (Khh) - ( B hhh) 2 « (l + S hl2 +5 hl3 +6 hh +26 hh S hl3 )C h C h C l3 (29) 

in the approximation of weak non-Gaussianity. 

The best estimates for the parameters are thus obtained by solving 

J2(B^,B^)f^ = (B^,B obs ) , (30) 
3 

1 We use the standard notation: (hhh) = [hhh + 5perms]/3!. 



while the statistical error on the parameters is deduced from the second-order derivatives of x 2 , 
which define the Fisher matrix, given here by 



Fij = {B® t B®). (31) 

For a real experiment, and if E-polarization is included as well, the above equations remain 
valid, except that the definition of the scalar product has to be replaced by a more complicated 
expression (see [5] for details). 

For each of the four isocurvature modes, the corresponding Fisher matrix has been computed 
numerically, including the polarization, for the noise characteristics of the Planck satellite in 
[U [5] . The error on the parameters f can then be deduced from the components of the Fisher 



matrix, via the expression 

For the various cases, we have obtained 

Af = {9.6,7.1,160,150,180,140} (CDM isocurvature) (33) 

Af = {9.6, 35, 4000, 720, 4300, 16600} (baryon isocurvature) (34) 

Af = {28, 36, 190, 150, 240, 320} (neutrino density isocurvature) (35) 

Af = {25,22,85,81,77,71} (neutrino velocity isocurvature) (36) 



As one can see, the error on the first two parameters is much smaller than the last four parameters 
in the CDM isocurvature case. An explanation for this result is given in [5J. 



5. Illustrative example 

To illustrate the previous results, which are model-independent, it is instructive to consider a 
simple class of models based on the presence of a spectator light scalar field during inflation, 
dubbed curvaton. This curvaton acquires nearly scale-invariant fluctuations during inflation 
and, later, behaves as a pressureless fluid when it oscillates at the bottom of its potential, before 
decaying. 

Here, we allow the curvaton a to decay into both radiation and CDM with the respective 
branching ratios j r and "y c . Since, in general, CDM can already be present before the decay, 
we define the fraction of CDM created by the decay as f c = 7 c f2 (T /(fi c + 'y c {l a ), where the fi's 
represent the relative abundances just before the decay. 

As shown in [2], the "primordial" adiabatic and isocurvature perturbations, i.e. defined after 



the curvaton decay, can be written in the form (16), with 



N< = — N< =— (37) 

N s a =-{f c -r), N^ c = -^j [/ c (l — 2/c) — r] , (38) 
c* erf 



where r = 37 r J7o-/[(4 — fJ CT )(l — (1 — j r )£l a )] is assumed to be small, since significant non- 
Gaussianities arise only if r -C 1. 

Let us first discuss linear perturbations. It is useful to introduce the curvaton contribution 
to the total adiabatic power spectrum 3 = (N^) 2 /[(N^) 2 + (N§) 2 ], where = H/<j) is 
associated with the inflaton fluctuation, and ./V^ = 0. 3 is directly related to the correlation 
C = P^g/y/PgP^ = vS sgn(/ c — r). The isocurvature-adiabatic ratio, given by 

a =PS- (^) 2 - 9 ( 1 fc) 2 ^ (39) 
~ (Nj) 2 + (N<) 2 ~ T r) - (39) 
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Figure 2. Plots of the coefficients (left panel) and of g^i^ L (right panel) as functions of f c . 

I, J, K, L are specified in the figure for each line. Other parameters are fixed as £ = 1, A = 10~ 3 
and r = 10~ 5 . 



is constrained by CMB observations [9] to be small which requires either f c ~ r or H <C 1 



Let us now turn to non-Gaussianities. Using (37), one finds /££ C = 3H 2 /(2r). This is the 



dominant contribution in the regime f c ~ r, the other components being suppressed. We thus 
concentrate on the more interesting case H <C 1 to discuss the size of the various components in 
terms of f c and r, considered as free parameters in our phenomenological approach. 

In the regime / c < r C 1, the purely adiabatic coefficient is the smallest one. The other 
ones are enhanced by powers of (—3) (since Nj~ /N^ = N^ a /N^ a = —3): 
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(40) 



where p is the number of U S" in the triplet {I,J,K}. In particular, the purely isocurvature 
coefficient is enhanced by a factor 27, but with the opposite sign: f^f S = —a 2 /(2r). All 
coefficients can be significant if r is sufficiently smaller than a 2 . 

In the opposite regime r < / c < 1, the purely adiabatic coefficient is, once again, the smallest 
one. All the coefficients are now positive and enhanced by factors (3/ c /r) p , where p is again the 
number of "5" indices: 



f ,JK _ (3f c V %CC K,CC _ « 2 r 3 
Jnl ~ \ r J /nl ' /nl ~ 54/4 



(41) 



Note that the enhancement factor is much bigger than in the previous case (40). The purely 
isocurvature coefficient, which dominates, is = a 2 / (2/ c ) and can be large if f c is sufficiently 

small, while the relative size of the other coefficients depends on the ratio r/f c . The full 
dependence of the coefficients on the parameter f c is illustrated in the left panel of Fig. [2j 

A similar analysis for the trispectrum has been presented in [3j. The usual coefficients tnl 
and #nl that describe the local trispectrum are then generalized into, respectively, nine TNL-like 
coefficients and eight ^NL-hke coefficients: 
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(42) 



The behaviour of the <?NL-like coefficients as functions of f c is plotted on the right panel of Fig.[2j 
The hierarchies are very similar to those observed for the bispectrum parameters. 

In conclusion, the above results show that a small isocurvature fraction in the power 
spectrum is compatible with a dominantly isocurvature non-Gaussianity detectable by Planck 
(e.g. a ~ and r <C f c — 10~ 8 yields f^f S — 5 x 10 3 ). Of course, the relations (40) or 

( |41[ ), are specific to the models considered here and would be a priori different in other models. 
It is therefore important to try to measure these six coefficients separately, in order to obtain 
model-independent constraints from observations. 
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